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In an earlier paper [K. K. Saha and A. Mookerjee, Phys. Rev. B 70 (2004) (in press) or, 
cond-mat/0403456 we had presented a formulation for the calculation of the configuration-averaged 
optical conductivity in random alloys. Our formulation is based on the augmented-space theorem 
introduced by one of us [A. Mookerjee, J. Phys. C: Solid State Phys. 6, 1340 (1973)]. In this 
communication we shall combine our formulation with the tight-binding linear muffin-tin orbitals 
(TB-LMTO) technique to study the optical conductivities of two alloys CusoAuso and NisoPtso. 

PACS numbers: 71.23.-k 



I. INTRODUCTION 

In this communication we propose to study the optical 
properties of disordered CuAu and NiPt 50-50 alloys from 
a first principle approach. We have chosen these two alloy 
systems because of several reasons : for CuAu, the bunch 
of d-like states sits about 1 eV below the Fermi level. 
For low photon energies, therefore, optical conductivity 
is dominated by the intra-band transitions within s — p- 
like states, which are extended and rather free electron- 
like. As a consequence, the optical conductivity for low 
photon energies below ~ 1 eV should have a Drude like 
behaviour. For higher photon energies inter-band transi- 
tions between the occupied d-states and the higher unoc- 
cupied states begin to take over. In sharp contrast, the 
Fermi energy of NiPt almost straddles the d-like peak. 
For this alloy the Drude behaviour should be confined to 
a very narrow low photon energy range. This contrasting 
behaviour should be reflected in our results. Moreover, in 
both the two alloy systems there is a large size mismatch 
between the constituents. This indicates that the stan- 
dard single site mean-field theories would be inadequate 
to capture the effect of this large size-mismatch. These 
alloy systems are therefore ideal to illustrate the advan- 
tages of the augmented-space recursion (ASR) method 
proposed by usi. 

Earlier theoretical work on optical conductivity for ran- 
dom alloys began with VelickjiS based on the single site 
coherent potential approximation (CPA) in an empiri- 
cal tight-binding model alloy. Butler— extended the ideas 
and combined the CPA with the first-principles Korringa- 
Kohn-Rostocker (KKR) technique. Banhart 4 used the 
KKR-CPA to study the optical conductivity of AgAu al- 
loys. This alloy system has close resemblance to CuAu. 
Banhart found discrepancies of his theoretical results 
with experiment^ and argued that various factors could 
be responsible : use of the density functional and the 
single-site mean-field approximations in theory and ef- 
fects of surfaces, their roughness, possible adsorbates and 
presence of large stresses in the samples, in the experi- 
ments. There have been a few more theoretical studies 
of optical properties of random alloys : Rhee et ali on 



CoAl, Uba et al. 8 on CoPt and Rhee et al& on Ni 3 Al. 
These works all base their approach on a large super-cell 
method to take care of the disorder. The method is brute 
force and less satisfactory than the CPA or ASR. 

In an earlier paper— we had presented a formulation for 
obtaining the configuration averaged optical conductiv- 
ity for random binary alloys. We had chosen as our basis 
the minimal set of the tight-binding linear muffin-tin or- 
bitals methocUi* 1 ^. Configuration averaging over various 
random atomic arrangements had been carried out us- 
ing the augmented-space formalism (ASF) introduced by 
us earlier for the study of electronic properties of dis- 
ordered system a 13 : 14 . The ASF goes beyond the usual 
single site mean-field approaches and takes into account 
configuration fluctuations about the mean-field. We shall 
present here a summary of the results derived in the ear- 
lier paper—. In linear response theory, at zero tempera- 
ture, the real part of the optical conductivity of a disor- 
dered alloy is given by the Kubo-Greenwood expression : 

a{u) = (1) 

where, the configuration averaged current-current corre- 
lation function <C S(ui) ^> is given by : 

i-^Tr| dE^lm{G^E)}i\ Im{G%E + w )}) (2) 
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If we define 

« S 1 (z 1 ,z 2 ) »= Tr ( j 7 G"( Zl ) jt G c (z 2 )) (3) 

then, using the Herglotz properties of the Green function, 
the correlation function becomes 

<S(w)»=t^-X! f dE [S 7 (E-,E+ +uj) 

71 7 ■* 

+ 5 7 [E+ , E~ + uj) - S y {E+ , E+ + u) 

-S 7 (E-,E-+lo)] (4) 

where 

/(£±) = limJ(E±iS). 
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The formulation we presented in our earlier paper— was 
based on the disordered induced scattering leading to the 
renormalization of the propagator G as well as the cur- 
rent terms j. The dominant contributions to the configu- 



J 



ration averaged correlation function was due to joint fluc- 
tuations of an averaged current and two electron propa- 
gators : 



f r/ 3 k 

« S^(z 1 ,z 2 ) » = / t-o < j 7 (k) » « G l '(k, z x ) » « j 7 (k) »t < G c (k, z 2 ) » 
Jbz ott 



(5) 
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We have also shown that disorder scattering renormalizes tor beyond the CPA approximation. We expressed the 

that averaged current to an effective current term and the effective current in terms of the self-energy as : 
average propagator to a configuration averaged propaga- 

I 



jf(k,Z!,Z 2 ) = «j 7 (k)»+2 



E(k,z 2 ) f(z 2 ) j 7 1)f (k) + j 7 1}t (k) f{ Zl ) £(M0 



+ E(k,za) f(z2) j 7 2) (k) f(zi) S(k,«i). 

and the contribution to the effective current to the correlation function was : 

d 3 k 



«SW(z l7 z 2 )» 



Tr 



jf (k, zx, 22) « G"(k, zi) » jf (k, zx, z 2 )t « G c (k, z 2 ) » 



'BZ 8^ 3 

The contribution of joint fluctuations between the two current terms and one propagator was given by 

d 3 k 



< S( 2) (z 1 ,z 2 ) > 



Tr 



BZ 



jW(k) f(zx) E(k,Z!) f(zx) jW^k) « G(k,z 2 ) » 



+ j« f (k) f(z 2 ) E(k,z 2 ) f(z 2 ) j«(k) « G(k,Z!) » 



The vertex correction terms in the ladder approximation contributed : 

«5 7 addor (z l! z 2 )» = Tr^T r 7 lL2 (^>^) f$* L *(z 1 ,z 2 )= Trr 7 (z 1 ,z 2 )®f 7 (z 1 ,z 2 ) A(z l7 z 2 ). 

L1L2 L3L4 



r 



(6) 



(J) 



(8) 



(9) 



where 



r d 3 v 

/ — G(k,z 2 ) Jf (k, z 1; z 2 ) G(k,zi) =r 7 (z 1; z 2 ) 
Jbz 8tH 



'BZ 

d 3 k' 



G(k',Z!) Jf (k',z 1 ,z 2 )tG(k',z 2 ) =f 7 (z 1; z 2 ) 



and 



/" d 3 k 

^z^Oi^) =y Bz ^3 GL 3 L 4 (k,zi) G L2Ll (k,z 2 ) 



= W l\ 6 l,l 2 Sl 3 l 4 



Wt=F L {z 2 ) 



l" 

+B l »(z 2 ) G RL » RL ,{z 2 )]] F L ,( Zl ). 



f l {z) = V^y d (%^) / « VAl » 

B L (z) = {y-x)d (^f£) I « 1/A L » 



Here x, y are the concentrations of the component atoms, 
C,A are the standard potential parameters of the TB- 
LMTO and d(f) — /a — /b- These super-matrices in 
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Lowest energy Vegard's Law 
Alloy lattice const lattice const 

(A) (A) 

CugoAuso 7.31 7.26 
NisoPtso 7.09 7.03 



TABLE I: Lowest energy and Vegard's Law lattice constants 
for CuAu and NiPt. 
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FIG. 1: The configuration averaged joint density of states 
and correlation function for CuAu (50-50) alloy shown as a 
function of the photon energy. 
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FIG. 2: Averaged optical conductivity, ao{u>) = S /u, and 
the density of states for a CuAu (50-50) alloy. The corrections 
to the optical conductivity are shown in Eqs. Ijril9|l . 

clear that the transition rate is dependent both on the 
initial and the final energies, throughout the frequency 
range of interest. That is : 

S(w) ± \T\ 2 J( W ) 

where 



{L} space are written as A and w. The full ladder vertex 
may now be written as 

A(zi,z 2 ) = y + w A w + w Aw Aw + ... 

= - A(*i,a&) ^r 1 (io) 

Including all sorts of corrections the averaged correlation 
function is then : 

< S 1 (z 1 ,z 2 ) >=< S[f ) {z ll z 2 ) > + < S < f ) {z 1 ,z 2 ) > 

+ «s; addor (zi,z 2 )» 

II. RESULTS AND DISCUSSION 

We have begun our study with the self-consistent TB- 
LMTO-ASR calculation on NiPt and CuAu 50-50 alloys. 
We have minimized the energy with respect to the vari- 
ation in the average lattice constant for both the alloys. 
The table shows the lowest energy lattice constants and 
compares them with the averaged or Vegard's law results. 
As expected, because of the large size difference between 
the constituents there is a "bowing" effect which is most 
prominent at the 50-50 alloys. The lowest energy lattice 
constant for both the alloys is greater than the Vegard's 
law predictions. 

Figure ^ shows the comparison between the scaled joint 
density of states (JDOS) and the averaged correlation 
function for a CuAu (50-50) alloy. From the figure it is 



J(w)= J dE J |^Tr(G c (k,£)G"(k,£ + c)). 

The Fig. ^ also shows that the disorder corrections to 
the current and the vertex correction are rather small 
and become negligible beyond photon energies of about 
22 eV. 

Figure [21 shows the optical conductivity <t(ui) for CuAu 
(50-50) alloy. The inset shows the configuration aver- 
aged density of states for the same alloy. The edge of 
the d-band complex is clearly seen to lie about 1 eV be- 
low the Fermi energy. The optical conductivity rapidly 
decreases as we increase the photon energy from zero up- 
wards. This decrease continues until about 1 eV and then 
the conductivity rises again and has considerable struc- 
ture as also shown in the correlation function for these 
photon energies (Fig.QJ. 

Figure shows the optical conductivity with a Drude fit 
[a D (u>) = <t(0)/(1 + {lot) 2 ), with o-(0) = 0.11, r = 9.78] 
for the lower photon energies. The Drude fit is good 
for photon energies below 1 eV. From this information 
we may deduce that for low photon energies the con- 
ductivity arises due to intra-band transition between the 
s — p states, which are free electron like and lead to a 
Drude type behaviour. Above 1 eV there is a onset of 
inter-band transition between the d and the conduction 
states and this leads to a sharp increase of optical con- 
ductivity and structure reflecting the structures in the 
d-like states. The inter-band contribution to the imagi- 
nary part of the dielectric function £2(0;) may be obtained 
from the optical conductivity data, by subtracting away 
the Drude contribution and dividing the result by u> : 
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FIG. 3: Averaged optical conductivity showing a Drude fit 
at low photon energies. 
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FIG. 4: Inter-band contribution to the imaginary part of the 
dielectric function for CuAu (50-50) alloy. 



£2(0;) = (rx(cj) — a D (ui))/ui . Below the onset of the inter- 
band transitions, this quantity vanishes. It then reaches 
a maximum at around 3 eV before decreasing. We have 
experimental data on AgAu (50-50)£, whose density of 
states closely resembles CuAu. The experimental data 
are in good qualitative agreement with Fig. The gen- 
eral shape with a shoulder around leV, a maximum and 
around 3 eV is clearly reproduced. 

Figure |S] shows the joint density of states and the av- 
eraged correlation function for the NiPt (50-50) alloy. 
The energy-frequency dependence of the effective transi- 
tion rate is considerable more pronounced than for CuAu. 
Disorder correction to the current terms and vertex cor- 
rections are also more in the low photon energy region. 
They become negligible for high photon energies. 

Figure El shows the density of states and the averaged 
optical conductivity for NiPt. Although the density of 
states for NiPt qualitatively resembles that for CuAu, 
unlike the latter, the Fermi level sits right atop the high 
peak due to the d-like states. The inter-band transitions 
between the ri-states and the conduction band is expected 
to start for very small photon energies, with a Drude 
contribution confined to a very narrow energy range near 
zero. The optical conductivity falls sharply in a very 



FIG. 5: The configuration averaged joint density of states 
and correlation function for NiPt (50-50) alloy shown as a 
function of photon energy. 
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FIG. 6: Averaged optical conductivity and the density of 
states for NiPt (50-50) alloy. 



narrow energy range and recovers almost immediately. 
This is expected from the density of states picture. Since 
the Drude ht is in a very narrow range indeed we do not 
show it explicitly in the figure. 

In Fig. we show the inter-band contribution to imagi- 
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FIG. 7: Inter-band contribution to the imaginary part of the 
dielectric function for NiPt (50-50) alloy. 
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nary part of the dielectric function for NiPt (50-50) alloy. 
The inter-band contribution begins at a very low photon- 
energy as expected and attains a maximum around 1 
eV. This is in contrast to the behaviour of CuAu, where 
Drude behaviour persists over a longer energy interval. 
We were unable to locate experimental data for this alloy 
system for comparison. 



III. CONCLUSION 

The frequency dependent transport quantities of two dis- 
ordered metallic alloys have been investigated applying a 
TB-LMTO-ASR based first principle theory 10 . We have 
found that the energy-frequency dependence to the effec- 
tive transition rate is more pronounced for NiPt than that 



of CuAu. Same tendency was observed in disordered cor- 
rection to the current terms and vertex corrections in the 
low photon energy region. We have also found that the 
conductivity occurs because of both intra-band and inter- 
band transitions. So the imaginary part of the dielectric 
function calculated from the conductivity splits into intra 
and inter-band contributions. For CuAu the intra-band 
transition takes place up to ~ 1 eV photon energy and 
in this region the conductivity curve follows Drude law 
of free electron model. For NiPt alloy the Drude be- 
haviour is confined to a very narrow energy range, as 
the Fermi energy straddles the d-like peak in density of 
states. We have compared our results with available the- 
oretical and experimental results and achieved a very im- 
pressive agreement with them. 



* kamal@bose.res.in 
t abhijit@bose.res.in 

1 T. Saha, I. Dasgupta and A. Mookerjee, J. Phys.: Condens. 
Matter 8 1979 (1996). 

2 B. Velicky, Phys. Rev. 184 614 (1963). 

3 W. H. Butler, Phys. Rev. B 31 3260 (1985). 

4 J. Banhart, Phys. Rev. Lett. 82 2139 (1999). 

5 J. Rivory, Phys. Rev. B 15 3119 (1977). 

6 P. O. Nielsson, Phys. Kondens. Matter 11 1 (1970). 

7 J. Y. Rhee, Y. V. Kudryavtsev, K. W. Kim and Y. P. Lee, 
J. Appl. Phys. 87 5887 (2000). 

8 L. Uba, S. Uba, V. N. Antonov, A. N. Yaresko and R. 



Gontarz, Phys. Rev. B 64 125105 (2001). 
9 J. Y. Rhee, Y. V. Kudryavtsev and Y. P. Lee, Phys. Rev. 
B 68 045104 (2003). 

10 K. K. Sa ha and A. Mookerjee , Phys. Rev. B 70 (2004) (in 
press) or, cond- mat/0403456 

11 O. K. Andersen, Phys. Rev. B 12 3060 (1975). 

12 O. Jepsen and O. K. Andersen, Solid St. Commun. 9 1763 
(1971). 

13 A. Mookerjee, J. Phys. C: Solid State Phys. 6 1340 (1973). 

14 S. Ghosh, N. Das and A. Mookerjee, Int. J. Mod. Phys. B 
21, 723 (1999). 



